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. Abstract. We show how Viennot's combinatorial theory of orthogonal polynomials 

^SJ ' may be used to generalize some recent results of Sukumar and Hodges on the matrix 

^ , entries in powers of certain operators in a representation of su(l,l). Our results 

Q ' link these calculations to finding the moments and inverse polynomial coefficients of 

, certain Laguerre polynomials and Meixner polynomials of the second kind. As an 

immediate consequence of results by Koelink, Groenevelt and Van Dcr Jeugt, for the 
related operators, substitutions into essentially the same Laguerre polynomials and 
Meixner polynomials of the second kind may be used to express their eigenvectors. 
Our combinatorial approach explains and generalizes this "coincidence" . 

< 

^ ■ Introduction 

-4— » 

In two recent papers (Hodges & Sukumar 2007; Sukumar & Hodges 2007), Sukumar 
and Hodges introduced a one-parameter family of operator algebras exhibiting a parity- 
^ , dependent structure, to study a quantum harmonic oscillator. The operators R and L, 
CN \ arising as the half of the square of the creation and annihilation operators, respectively, 
generate a subalgebra representing the Lie algebra sm(1,1). This representation is 

■ the direct sum of two discrete series representations. Sukumar and Hodges observed 
! that calculating the matrix entries in the powers oi L + R lead to some interesting 

O \ combinatorial questions. In this paper we show that these results may be generalized 
in such a way that the answer may be expressed in terms of moments and inverse 
polynomial coefficients associated to Laguerre polynomials and Meixner polynomials of 
the second kind. 

■ Almost at the same time, Klimyk (2006) used discrete series representations of 
sm(1, 1) to study another model of a quantum harmonic oscillator. In Klimyk's set- 
ting, the operator L + R corresponds to the momentum operator, and Klimyk (2006) 
noted that the eigenfunctions of this operator may be expressed via substitutions into 
Meixner-Pollaczek polynomials with the appropriate parameters. This observation is 
a consequence of the more general results of Van Der Jeugt (1997), Koelink & Van 
Der Jeugt (1998) and Groenevelt & Koelink (2002), describing the spectra of certain 
self-adjoint elements in an arbitrary discrete series representation of su(l,l), where 
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besides the Meixner-PoUaczek polynomials, the Laguerre polynomials and the Meixner 
polynomials of the first kind also appear, depending on the choice of the parameters. 

In this paper we generalize and explain this remarkable "coincidence", using Vien- 
not's (1983) combinatorial theory of orthogonal polynomials. For that purpose we will 
consider tridiagonal operators T defined on a subspace of P{Z>q) that contains at least 
the subspace generated by the basis vectors {cq, ei, 62, . . .}. 

In section E] we express the matrix entries (ej+rf,T™ei) in terms of weighted Motzkin 
paths. Using Viennot's (1983) theory, in the case when the matrix of T has nonzero 
off-diagonal entries wherever this is allowed, we are able to associate an orthogonal 
polynomial sequence {pn{x)}n>o to T in such a way that the matrix entries (eo,T™'eo) 
and (crf, T'^cq) may be expressed in terms of the moments and inverse polynomial co- 
efficients associated to {p„(a;)}„>o. Applications of these general results to the discrete 
series representations of su{l, 1) follow in section [31 

In section H] we see that for those tridiagonal operators T that have nonzero off- 
diagonal matrix entries wherever this is allowed, the same orthogonal polynomial se- 
quence {pn{.x)}n>o as in section [2] may be associated to T to express its (potential) 
eigenvectors. For the discrete series representations of su(l, 1) we thus recover the ne- 
cessity part of the above mentioned results of Van Der Jeugt (1997), Koelink & Van 
Der Jeugt (1998) and Groenevelt & Koelink (2002). 

The ideas presented in this paper also provide a combinatorial framework to discuss 
some representations of other Lie (super) algebras over £^(Z>o). In the concluding sec- 
tion [5] we outline a sample application. The other example in section [5] indicates that 
we do not have to limit ourself to self-adjoint operators at this level of investigation. 
It seems harder to replace the Hilbert space £^(Z>o) with another Hilbert space, as in 
that case one would be facing the dilemma, what notion should replace the notion of 
the associated orthogonal polynomial sequences. Finding such a generalization seems 
an interesting challenge for future research. 

1. Preliminaries 

1.1. Viennot's combinatorial proof of Favard's theorem. Concerning orthogonal 
polynomials, in this paper we follow Viennot's (1983) notation and terminology, but 
sometimes we complement the facts stated therein with results cited from Chihara's 
(1978) classical work. Viennot's work (1983) is currently out of print but available on 
the author's webpage. Some of his most important results were also communicated by 
Stanton (2003), a recommended source for those readers whose French is not fluent. 

The direct way to define an orthogonal polynomial sequence (OPS) {pn{,x)}n>o is to 
provide a linear form / : C[a;] — > C and postulate the following three axioms: 

(i) for all n, p„(x) is a polynomial of degree n. 
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(ii) f {pm{x)pn{x)) = if m 7^ n, 

(iii) for all n, f{pn{xf) 7^ 0. 

The map / is called a moment functional. Whenever an OPS exists, each of its elements 
is determined up to a non-zero constant factor, see the corollary of theorem 2.2 in ch. 
I of Chihara's (1978) book. 

An equivalent way define an OPS is by means of Favard's theorem, see theorems 
4.1 and 4.4 in ch. I of Chihara's (1978) book. This states that a sequence of monic 
polynomials {pn{x)}n>o is an orthogonal polynomial sequence, if and only if it satisfies 
the initial conditions 

Po{x) = l,pi{x) = X -bo, (1.1) 
and a two-term recurrence formula 

Pn+l{x) = (X - bn)Pnix) - A„p„,_i(x) fom > 1, (1.2) 

where the numbers 6„ and A„ are constants and A„ 7^ for n > 1. The above form 
appears in theorem 9 of ch. I in Viennot's (1983) book, in Chihara's (1978) work the 
indices are shifted. The original proof of Favard's theorem provides only a recursive 
description of the moment associated functional /. In proposition 17 of ch. I, Viennot 
(1983) gave an explicit combinatorial description by expressing the moments /i„ := 
/(x") as sums of weighted Motzkin paths. A Motzkin path of length n, as defined in 
def. 15 of ch. I by Viennot (1983), is a path oj = (sq, . . . , s„) in Z x Z from sq to s„ 
such that the second coordinate of all Sj's is non-negative, and each step (sj,Sj+i) is 
either a northeast step (1, 1) or an east step (1, 0) or a southeast step (1, —1). Viennot 
(1983) introduces the following valuation: the weight v{uj) of uj is the product of the 
weights of its steps (sj_i, Si), where each northeast step has weight 1, each east step at 
level k has weight bk and each southeast step starting at level k has weight A^. Here 
the level is the second coordinate of the lattice point. He then defines /i„ := J2u)'^i^) 
as the total weight of all Motzkin paths uj from (0, 0) to (n, 0). 

Theorem 1.1 (Viennot (1983)). Let {pn{x)}n>o be a sequence of monic polynomials, 
given by U.l\) and U.^) and let f : C[x] C be the linear map given by f{x^) := iin- 
Then for all n,k,i > we have 

f{x'^Pk{x)pi{x)) = Ai • • ■ A^ ^ v{uj) 
where the summation is over all Motzkin paths of length n from level k to level i. 

Substituting n = in theorem 11.11 yields the more difficult implication of Favard's 
theorem. A generalization of theorem 11.11 is theorem 1 in ch. Ill of Viennot's (1983) 
work, which allows the computation of the inverse polynomials of a system of monic 
polynomials {p„(x)}„>o given by f 1 1.11) and fll.2p . even if some scalars A„ are zero (and 
thus {pn{x)}n>o IS uot an OPS). The inverse polynomials g„(x) := YTi= =0 1n,iX^ are 
defined by = Y^^^^o 9n,iPi{x) . 

Theorem 1.2 (Viennot (1983)). Let {pn{x)}n>o be a system of monic polynomials 
defined by U.l\) and U.^) for some numbers {6„}„>o and {A„}„>i. The coefficient qn,k 
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0/ in the inverse polynomial qn{x) is then the total weight of all weighted Motzkin 
paths of length n starting at (0, 0) and ending at {n, k). 

1.2. Viennot's "histoires" , Laguerre and Meixner polynomials. Sometimes mod- 
els of even deeper combinatorial interest may be built, using Viennot's (1983) second 
valuation of Motzkin paths, defined as follows. Given the sequences of numbers {ciA;}fc>0) 
{bk}k>o and {ck}k>i, we define the weight Vi{uj) of a Motzkin path u as the product of 
the weights f i(si_i, Si) of its steps, such that each northeast (resp. east, resp. southeast) 
step starting at level k has weight (resp. b^, resp. c^). Setting 

Afc = ctfc-iCfc for A; > 1, (1.3) 

for a Motzkin path u from (0, 0) to (n, 0), Vi{uj) is equal to v{u), as defined in subsec- 
tion 11.11 because each northeast step starting at some level k — 1 may be matched to 
the first subsequent southeast step starting at level k. More generally, for a Motzkin 
path uj starting at level k and ending at level /, we have 

{E<^^iH if k = l, 

akUk+i- ■ ■ai^iY.c.M^) ifk<l, (1.4) 
CkCi-i- ■ -ci+iJ^ujM^) iik>l, 

see lemma 1 in ch. II of Viennot (1983). Algebraically, eq. (11. 3p corresponds to switch 
from the three term recurrence relation for monic orthogonal polynomials to a gen- 
eral three term recurrence relation. The combinatorial interest in vi arises when 
the sequences {ak}k>o, {bk}k>o and {ck}k>i consist of non-negative integers, allow- 
ing us to think of these weights as making choices from a set of options. In such 
situations Viennot (1983) defines a story ("histoire") as a pair (cu; (pi, . . . ,p„)) of a 
Motzkin path u = (sq, . . . , s^) and a sequence {pi, . . . ,Pn) of positive integers satis- 
fying 1 < Pi < vi{si-i,Si) for 1 < i < n. Clearly the moment /i„ is the number of 
stories of length n. Thus, in such situations, we may replace the weighted lattice path 
model with a model that involves enumerating (non-weighted) combinatorial objects. 
In particular, permutation enumeration models of the (monic) Laguerre polynomials 
and of the Meixner polynomials of the second kind (discussed below) may be obtained 
by specializing Viennot's (1983) bijection between his "histoires de Laguerre" of length 
n associated to the valuation 

ak = k + 1, bk = 2k + 2 for > 0; Ck = k + 1 for > 1, (1.5) 

and the permutations of the set {1, . . . ,n + 1}. For the details we refer to Viennot's 
(1983) work. 

The (monic) Laguerre polynomials Li^\x) are the OPS defined by (11. ip and (II. 2p 
where bn = 2n + a + 1 and = n{n + a), see §5 in ch. II of Viennot's (1983) work. 
The associated moments are given in (31') of ch. II by Viennot's (1983): 

/i„ = (a + 1)„ = {a + l){a + 2) ■ ■ ■ (a + n). (1.6) 

The combinatorial model mentioned above is associated to the case a = 1 yielding 
large Laguerre stories. For general a we have weighted Laguerre stories, for a = we 
obtain restricted Laguerre stories by limiting 6'^ to k above. The restricted Laguerre 
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stories are thus a subset of the large Laguerre stories and correspond bijectively to the 
permutations a G Sn+\ with (t(1) = n + 1. 

Another class of particular interest to us are the Meixner polynomials of the first kind 
mn{x;6,ri) := (c — l)"/c"m„(x; /5, c), whose monic variant {m„(a;; /5, c)}„>o is defined 
by (mH) and where 

(1 + c)n + (3c cnjn + (3-1) 

On = and Xn = -. . (1.7) 

1 — c (1 — c)^ 

Viennot (1983) (see (54") in ch. 2) )has a combinatorial proof of the fact that the 
moments of the corresponding linear functional are given by 

MAc) = (l-c)^5^A:V^. (1.8) 

fc>0 

Here {p)k = P{P + I) ■ ■ ■ {P + k - 1). 

Finally we will be interested in the (monic) Meixner polynomials of the second kind 
Mn{x; S,ri), defined by (11. ip and (11.21) where 

bn = i2n + ri)5 and A„ = (5^ + l)n(n + r/ - 1). (1.9) 

A combinatorially interesting special case is 5 = 0, implying 6„ = and A„ = n{n+r]—l). 
The moments associated to {M„(a;; 6, ri)}n>o may be expressed using the Motzkin paths 
associated to the Laguerre polynomials {Ln~^\x)}n>o subject to the restriction that 
no east step occurs. Motzkin paths with no east steps are also known as Dyck paths. 
Using the permutation-enumeration model mentioned above, Viennot (1983) (in ch. 2, 
example 21) shows that fin{0,ri) is given by 

/in(0,r^)= 5^r^^('^). (1.10) 

Here Zn is the set of alternating or zig-zag permutations of {1, 2, . . . , n, }, starting and 
ending with a rise, and s(cr) is the number of left-to-right minima ("elements saillants") 
a(i), defined by (t(z) = min{(T(l), cr(2), . . . , cr(2)}. A generating function for the numbers 
fJ'2n{0,v) for positive integer rj was already computed by Carlitz (1959). He considered 
the polynomials fn\x) = (— i)"M„(ix; 0, A;) for i = ^/^, but using the substitution 
X H- >■ x/i we may recover the moments /U2„(0,?7). In the next lemma we state Carlitz' 
(1959) result (see his formulas (9.3) and (9.13)) modified for our purposes, together 
with the outline of a combinatorial proof that extends its validity to all numbers rj. 



Lemma 1.3. For all rj ^ the moments of {Mn{x;0,ri)}n>o satisfy 
f,2n{0,v) = Ei'J where sec''(t) = ^ — -. 

n>0 ^ 



Proof. Let E2n,k denote the number of those zig-zag permutations on 2n elements that 
start and end with a rise and have k left-to-right minima. Since 1 is the rightmost 
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left-to-right minimum, using the decomposition cr(l) ■ ■ ■ a{2n) = cr(l) ■ ■ • 1 ■ ■ ■ a{2n) we 
have the recursion formula 




(1.11) 



Here T2n-2m-i counts the number of zig-zag permutations on 2n — 2m — 1 elements, 
starting with a descent and ending with a rise. Introducing 0fc(t) := X]n>o ^2n,kt'^"'/ (2n)! 
we obtain the recursion ^ 

4>k{t) = / (f)k-i{u)tan{u) du. 
Jo 

Using this, we may show by induction that (pk(t) = (— lncos(t))^/A;!. Thus, by (11.101) . 
f^2n{0,ri) is the coefficient of t^"/(2n)! in 

^(-l)''lncos(t)'' • 77''//;;! = exp(-?7lncos(t)) = sec(t)''. 

A;>0 

□ 



A variant of the Meixner polynomials of the second kind are the Meixner-Pollaczek 
polynomials {Pn^\x; (j))}n>o, defined by (1.7.1) in the Askey- Wilson scheme by Koekoek 
& Swarttouw (1998). We should note for future reference that the two OPS are essen- 
tially the same. The formula connecting them is (3.22) in ch. VI of Chihara's (1978) 
work: 

sin^ (b 

P^{x;(p) = —^Mn{2x;6,2\), 5 = cot0, < < vr. (1.12) 

Chihara (1978) does not use the adjective "Meixner-Pollaczek", but gives the same re- 
currence for them in (5.13) of his ch. VI as the recurrence (1.7.3) of Koekoek & Swart- 
touw (1998). Setting = n/2 in (fTT^ leads to 5 = 0. According to (1.7.4) of Koekoek 
& Swarttouw (1998), the normalized version of P^{x](f)) is n!/(2 sin0)"P^(a;, 0). As a 
consequence of eq. f ll.l2p we obtain the following. 

Corollary 1.4. For < < 1, the normalized version of the Meixner-Pollaczek poly- 
nomials {P^(x;0)}„>o is {2~"M„(2x, cot 0, 2A)}„>o. 



Finally we review the inverse polynomials of {Ll^\x)}n>o, {rhn{x; (3,c)}n>o and 
{Mn{x; 6,ri)}n>o. All three OPS are examples of Sheffer orthogonal polynomials de- 
fined as polynomials {pn{x)}n>o, given by a generating function 

Vpn(x)^ = fit) expixgit)). (1.13) 

n>0 

As noted in (11) and (12) of ch. Ill in Viennot's (1983) work, the inverse polynomials 
{g„(x)}n>o of a Sheffer OPS {pn{x)}n>o form a Sheffer OPS, given by 

j.n 1 



n>0 



Here g^ ^^(t) stands for the compositional inverse of g{t). 
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1.3. Discrete series representations of su{l, 1) and two applications. In this 
subsection we review the discrete series representations of the Lie algebra 
together with two applications of them in models of a quantum oscillator. Our notation 
will be close to the one used by Groenevelt & Koelink (2002). 

The Lie algebra su{l, 1) has the generators H, B, C, subject to the commutation 
relations [H, B] = 2B, [H, C] = -2C and [B, C] = H. According to section 6.4 in 
the work of Vilenkin and Klimyk (1991), there are four classes of irreducible unitary 
representations of sm(1, 1): the positive and negative discrete series representations, the 
principal series representations and the complementary series representations. In this 
paper we will focus on the two discrete series representations, defined on the represen- 
tation space £^(Z>o), whose orthonormal basis vectors we will denote by {e„}„>o. The 
positive discrete series representations vr^ are labeled by A; > 0. The action is given by 

■K^{H)en = 2{ k + n)e„, 

7r+(5)e„ = ^/{ n + l){2kTn}e n+i and (1.15) 
7r+(C)e„ = -^/n{2kTr^^en-l, 

see eq. (2.2) of Groenevelt & Koelink (2002). The negative discrete series representa- 
tions TT^ are labeled hj k > 0. The action is given by 

7r^(if)e„ = -2{ k + ra)e^, 

T^~^{B)en = —^/n{2k + n — l)e„_i and (1-16) 
7r^(C)e„ = ^/inTT)i2kTn^en+l, 

see eq. (2.3) of Groenevelt & Koelink (2002). The principal series and complementary 
series representations are defined on £^(Z), a different Hilbert space. Koelink & Van 
Der Jeugt (1998) (proposition 3.1, rephrased) have expressed the eigenvectors of vr^(c- 
H + B — C) for an arbitrary c G R satisfying |c| < 1 in terms of the orthonormal 
Meixner-Pollaczek polynomials p^\x; 0), given by 

where the polynomials P^\x; (p) are the Meixner-Pollaczek polynomials. 

Theorem 1.5 (Koelink & Van Der Jeugt). Introducing Xfp = — cos(f)H + B — C , the 
vector 

n>0 

is a generalized eigenvector for n^i^X^) , with eigenvalue 2xsm.(j). 

An analogous result for c = 1 was worked out by Van Der Jeugt (1997), leading to 
Laguerre polynomials, who also mentions (in his section VIII) that the case |c| > 1 leads 
to Meixner polynomials of the first kind. The extension of Theorem 11.51 to negative 
discrete series representations is stated in proposition 3.1 of Groenevelt & Koelink 
(2002). An excellent source treating all cases simultaneously is ch. 3 of Groenevelt's 
thesis (2004). 
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Klimyk (2006) used the positive series representations of sm(1, 1) to study a model 
of a quantum oscillator. He defines 1) as the Lie algebra generated by Jq, J+, J-, 
subject to the commutation relations [Jo, J+] = J+, [Jo, J-] = — J- and [J_, J+] = 2 Jo. 
This definition is equivalent to the above one, via setting H := 2Jo, B := J+ and 
C := — J_. Klimyk (2006) realizes the orthonormal basis vectors {e„}„>o as polynomials 

1 /2 

^n{y) •= (■^^^^IT^) y" single variable y, and the Hilbert space £^(Z>o) is 

realized as the closure of the space of all polynomials in a single variable y, whereas 
the generators Jo, J+, J- of sm(1, 1) become differential operators. One of Klimyk's 

(2006) observations (eq. (17)) is that the eigenfunctions of the momentum operator 
Ji := ^(J+ + J-) are of the form 

My) = J2Pn\p;^my''. (1.17) 

n>0 

with the Meixner-PoUaczek polynomials Pn\p]T^/'^)- Obviously, fll.l7p may be ob- 
tained from Theorem 11.51 by substituting = 7r/2. By eq. (1.7.11) of Koekoek & 
Swarttouw (1998), eq. ffTTTD is equivalent to Vp(y) = (1 - iyY^^''^\l + lyY^''-^. 

Sukumar & Hodges (2007) considered another model of a quantum oscillator, using 
(besides others) the operators L, _R, S", acting on £^(Z>o) as follows: 

Le2n = v/ (2n + l + a)(2n + 2))/2 

Le2n-i = V(2n + l + a)(2n))/2 

Re2n+2 = V(2n + l + «)(2n + 2))/2-e2„ (1-18) 

i?e2„+i = v^(2n + 1 + «)(2n))/2-e2„-i 

S'e„ = (2n + 1 + a)/2 ■ e„ 

The parameter a G M is assumed to satisfy \a\ < 1. As noted by Sukumar & Hodges 

(2007) in eq. (2.1), the operators L,R,S obey the same commutation rules as the 
generators B, —C, H above, thus we obtain a representation tt of su{l, 1) by setting 

7t{B) := L, 7r(C) := -R, 7t{H) := S. (1.19) 

This representation may be written as a direct sum of two representations by considering 
the restrictions of the operators L, R, S to the closure of the subspaces generated by 
{e2n}n>o and {e2n+i}n>o, respectively. Introducing tt/j to denote the representation of 
su{l, 1) given by 

np{B)en = V( n + l)(n+^ 

vr/3(C)e„ = -^n{n-l + f3)en-i, (1-20) 

7rf3{H)en = 2{n + (3/2)en, 
it is not difficult to see that first restriction is equivalent to n(^a+i)/2, whereas the second 
restriction is equivalent to vr(Q+3)/2- Comparing fll.20p with fll.lSp we obtain vr^ = n'^^^ 
for f3>0. 

Corollary 1.6. For a G (—1, 1], the representation vr of su{l, 1) given by fll.l9p is iso- 
morphic to the direct sum of the positive discrete discrete series representations 'n'^a+i)/4 
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Remark 1.7. For /5 = 0, we get 7ro(-B)eo = 7ro(C)eo = 7ro(-ff)eo = 0, thus we may 
decompose ttq as a direct sum of zero acting on the hnear span of eo and of the restriction 
of 7Tq acting on the orthogonal complement Cq of Cq. By shifting basis vector indices in 
Cq we may show that the restriction of ttq to Cg is equivalent to 112 = n^. 



It should be emphasized that corollary 11.61 is applicable only to the representation of 
sm(1,1), induced by the operators L,R,S. In the work of Sukumar & Hodges (2007) 
these operators are expressed in terms of the annihilation operator A and its adjoint . 
The algebra generated by A and A'' is a representation of a Lie superalgebra properly 
containing sm(1, 1). In this paper we focus on the combinatorial statements related to 
the operators L, R, S and generalizations of these. In particular, Hodges & Sukumar 
(2007) show (see eq. (5.2)) that for a = 1 we have 

(eo, {L + Rf^eo) = E^m and (d, (L + Rf'^ei) = T^m+i for all m > 0. (1.21) 

Here the numbers {E2m}m>o resp. {T2m+i}m>o are the secant resp. tangent numbers, 
given by 

sec(x) = y and tan{x) = Y 

^ ' ^ 2m ! ^ ' ^ 2m + 1 ! 

m>0 ^ ' m>0 ^ ' 

At the light of corollary II. 6[ eq. fll.2ip is equivalent to 

(eo, 7r+ (5 - C)™eo) = and (eg, 7r+(5 - C)"eo) = Ts^+i. (1.22) 



2. Matrix entries of tridiagonal operators on £^(Z>o) 



Definition 2.1. Consider the complex Hilbert space £^(Z>o) with orthonormal basis 
vectors Cq, Ci, 62, .... We call a linear operator T : -D(T) C £^(Z>o) £^(Z>o), 
defined on at least the linear span of {eo, 61,62, ■■■ ,}, a tridiagonal operator, if there 
exists real numbers {ln}n>o, {dn}n>o and {un}n>o such that Tcq = doCo + IqCi and, for 
each n > 1, we have Tcn = M„e„_i + (i„e„ + /„e„+i. 

Equivalently, the restriction of T to the linear span of {eo, ei, 62, ■ ■ ■ ,} may be repre- 
sented by a tridiagonal matrix 

/ do uo 

Iq di ui 

li d2 U2 

I2 ds u-i 

\ : 

such that all entries are real numbers. 

The letters {ln}n>o, {un}n>o and {dn}n>o should remind the reader of the words 
"lower", "upper", and "diagonal", as it is customary in numerical analysis. We define 
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the LDU -decomposition of the tridiagonal operator T as the sum of operators T = 
L + D + U, where the operators L, D, U are given by f/eo = and 

Le„ = /„e„+i, [/e„+i = M„e„ De„ = rf„e„ forn > 0. (2.1) 

Remark 2.2. A it is a fortunate coincidence that the operator L used by Sukumar & 
Hodges (2007) corresponds to the operator L in this setting, modulo corollary II .61 The 
choice of D will depend on what we want to compute: to recover the formulas regarding 
the matrix entries of L + R and L + i^ + S" in the work of Sukumar & Hodges (2007) and 
Hodges & Sukumar (2007) we would need to equate D to zero or to the appropriate 
restriction of S. 

The following lemma may be considered as a generalization of equations (4.2) and 
(4.3) in Hodges & Sukumar (2007). A reader making this comparison should note that, 
for technical reasons, we read the words X consisting of the letters L, D, U right-to-left. 

Lemma 2.3. Let X he a word of length m consisting of the letters L, D, U, where 
T = L + D + U is the LDU -decomposition of a tridiagonal operator T : D{T) C 
£^(Z>o) £^(Z>o) and let i > and d> —i be integers. Associate to the pair (X, i) a 
lattice path starting at (0, i) such that each L is a southeast step, each U is a northeast 
step, each D is an east step, and we read the letters in X right-to- left. Then {ei^d,X^i) 
is not zero only if the lattice path associated to {X, i) is a Motzkin path of length m 
from (0, i) to (m, i -\- d). 

Proof. Applying X to Cj involves applying the operators L, D, U in it, reading X from 
right to left. While applying these operators, at every instance of the calculation our 
partial result is a scalar multiple of a single basis vector. Each L increases the index of 
this basis vector by 1, each U decreases it by 1 and each D leaves the index unchanged. 
Thus the end result is either zero or a scalar multiple of Ci+d where d is the difference 
between the number of Ls and Us in X (and may be negative). If at any instance 
the number of f/s read exceeds the number of Ls by more than i then we get zero 
hy Ucq = 0. Therefore, if {ei+d,X^i) then {X,i) must represent a Motzkin path 
starting at (0,i) and necessarily ending at (m,i + rf). □ 

The next lemma generalizes equations (4.6) and (4.7) of Hodges & Sukumar (2007). 
Our formulas look somewhat different mainly because of our convention of reading X 
right-to- left. 

Lemma 2.4. For allp,i > and d > —i we have 

{ei+d,XDLPei) = dp+i ■ {ci+d^XL^Ci) and (2.2) 
{ci+d.XUL'P-^^ei) = Ip+iUp+i ■ {ei+d,XLPei). (2.3) 

The proof is immediate. Using eq. (11.41) and lemmas [2.31 and [2l4l we may express each 
(cj+rf, T^Cj) in terms of a total weight of weighted Motzkin paths. 
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Theorem 2.5. Let i > and d > -i be integers, T : D{T) C e^{Z>o) ^^(Z>o) « 
tridiagonal operator and T = L + D + U its LDU -decomposition. Then we have 

[ui-iUi-2- ■ -Ud^^viuj) ifd<0. 

Here J2uj "^i^) total weight of all weighted Motzkin paths of length m starting at 

(0, i), ending at {m,i + d) such that each northeast step has weight 1, each southeast step 
starting at level n has weight A„ = /„_iM„_i, and each east step at level n has weight dn- 
In particular, if D is identically zero then Yluj'^i^) total weight of all weighted 

Dyck paths of length m starting at {0,i), ending at (m, i + d) such that each northeast 
step has weight 1 and each southeast step starting at level n has weight A„ = ln~iUn~i- 

As a consequence of theorem 11.11 setting d = Q and i = in theorem 12.51 yields the 
following algebraic statement. 

Corollary 2.6. IflnUn 7^ for all n> then for all m, (cq, T^^cq) is the mth moment 
of the functional associated to the OPS given by ( li.ij) and lil.B) where bn = dn and 

Remark 2.7. As noted in the Preliminaries, an OPS and the associated moments deter- 
mine each other only up to a constant factor. In corollary 12.61 we understand that the 
moments were defined as total weights of weighted Motzkin paths as in theorem II. 1[ 
If we obtain the moments by some other means, we always need to check whether an 
adjustment by a constant factor is necessary. It suffices to check whether /xq = 1 is 
satisfied. 

When we apply theorem 12.51 to c? 7^ 0, theorem 11.21 becomes useful, at least for the 
case i = 0. 

Theorem 2.8. For all m,d > we have 

where firn,d is the coefficient of x'^ in the degree m inverse polynomial for system of 
monic polynomials given by (QHP and ( (i.^) where bn = dn and A„ = In-iUn-i- 

3. Matrix entries in discrete series representations of su(1, 1) 

For a (positive, or negative) discrete series representation vr of sm(1, 1), the operator 

Te[7r,c] :=7r(-C) + c-7r(/J) + 7r(5) (3.1) 

is a tridiagonal operator for any c G M. Furthermore, the operators appearing in the 
[/-decomposition of T[7r, c] are 

L[7r] = 7r(-C), D[7r,c] = c- 7r(iJ), [/[tt] = 7r(fi). 
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It is easy to see that /„m„ 7^ holds in all cases, thus we may apply corollary 12.61 to 
compute (eo, (7r(— C) + c ■ n{H) + ■n{B))'^eo) for every discrete series representation tt, 
all c e M and all m G N. In particular, using fll.lSp and f ll.l6p we obtain following 
corollary. 

Corollary 3.1. For any A; > 0, m G N and c G M, the matrix entry (cq, (7r^(— C) + 
c ■ T[^{H) + 7rj!"(i?))'"eo) is the mth moment of the functional associated to the OPS 
c]„(x)}„,>o given by ( fi.ij) and ( fi.^j) where, 

bn = 2c{k + n) and Xn = n{n + 2k — 1). (3-2) 

Similarly, the matrix entry {cq, (tt^ (— C) + c ■ t:^{H) + {B))"^eo) is the mth moment 
of the functional associated to the OPS {p[k, — c]„(x)}„>o. 

The appropriate parts of the next proposition (restricted to the cases |c| = 1, |c| > 1 
and |c| < 1, respectively) are either implied or explicitly stated in some equivalent form 
in the works of Van Der Jeugt (1997), Koelink & Van Der Jeugt (1998) and Groenevelt 
& Koelink (2002). A detailed and deep simultaneous discussion of all three cases may 
be found in ch. 3 of Groenevelt's thesis (2004). Here we make a summary statement 
with respect to our bases and outline the proof that ensues from applying Meixner's 
classical method. 

Proposition 3.2 (Groenevelt-Koelink-Van Der Jeugt). The OPS {p[k,c\n{x)}n>o de- 
fined in corollaru \3 . 1\ is given by 



p[k,c]n{x) 



\v/^-c)"m„(^-fc;2fc,^^) ./|c|>l, 

lT-'\x) ^fc=l, 

irL^n''''\-x) ifc = -l, 

^{VT^)'' ■ Mn{x/VT^;c/VT^,2k) if\c\<l. 

- (2fc-l) 



Here the polynomials {Ln (a;)}n>o ore Laguerre polynomials, the polynomials {mn{x; 2k, (c— 
a/c^ — l)/(c + Vc^ — l))}n>o anc? {M„(x; c/\/l — c^, 2k)}n>o, respectively, are Meixner 
polynomials of the first and second kind, respectively. 

Proof. We apply Meixner's method to classify all OPS {Pn{x)}n>o given by (11.11) and 
where 

bn = dn + f and Xn = n{gn + h), (3-3) 

see section 3 of ch. 6 in Chihara (1978). By (13.21) we have 

d = 2c, f = 2ck, g = l, h = 2k-l. 

The case d'^ — ig = A{c^ — 1) > is thus equivalent to |c| > 1. Following section 3 of 
ch. 6 in Chihara (1978) we set p := \J — Ag = 2\J(? — 1. Note next that replacing 
c\n{x) with 

q\k, c]„(x) := c]„(x - 7^,^) 
yields an OPS {^ni^, c](x)}„>o satisfying (13. 3p with 

d = 2c, f = 2ck + 7fc,c, fi' = 1; h = 2k — 1. 
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We want to choose '^k,c in such a way that 



f = '^^j~!~^\ i.e., 2c/c + 7fc, - 



d + p ' ' ' 2c + 2V^^ 

is satisfied. This yields 

2k 



Ik,, 



2kc = 2k{c - Vc2 - 1) - 2kc = -2k^ - 1 



and we may apply the results stated in section 3 of ch. 6 in Chihara (1978) to c](a;)}„>o. 
The parameter "c" appearing in the definition of the corresponding Meixner polynomi- 
als of the first kind becomes 

d-p 2c -2^/1?^ c-^/W^ 



d+p 2c + 2V^^^ c+Vc^^ 
and the parameter (3 becomes 1 + (2k — 1)/1 = 2k. We obtain 
1 



or, equivalently, 



[k, c]n{2Vc^ - IX) = TJlnix; 2k, (c - Vc^ - 1)/(C + - 1)) 



q[k, c]„(x) = (Vc2 - 1 - c)"m„ ( J^— ] 2k, ^ 1 1 

The formula for p[k, c]„(x) now follows from p[k, c]„(x) = q[k, c]„(x + ■yk,c)- 

The case d"^ — Ag = 4(c^ ~ 1) < is equivalent to |c| < 1. The parameter S of 
the corresponding Meixner polynomials of the second kind becomes 2c/ V A- 4c2 = 
cj \J\ — and the parameter i] becomes 1 + {2k — 1)/1 = 2k. We obtain 

^ ' p[A;, c]„(xVl - c2) = M„(x; c/Vl - c^, 2k). 



The case d'^ — Ag = A{c^ — 1) = is equivalent to |c| = 1. For c = 1, (13. 2p gives 
bn = 2n + 2k and A„ = r;,(?7, + 2A; — 1), which define the monic Laguerre polynomials 
p[k, l]„(x) = Ln''~^\x). Similarly, we also have p[k, — l]„(x) = {—l)'^Ln''~^\—x). □ 

Remark 3.3. Inspired by theorem 11.51 and corollary 11.41 for |c| < 1 we may introduce a 
(p such that < < TT and c = cos(0) hold. Then we obtain 

p[k, cos(0)]„(x) = sin(0)''M„(x/ sin(0); cot(0), 2k). (3.4) 

As indicated in the preliminaries, a combinatorial interpretation of the "Laguerre 
case" |c| = 1 may be obtained using Viennot's (1983) "Laguerre stories". The "best" 
combinatorial interpretations (linked to permutation enumeration with no special weights) 
is associated to the case k G {0, 1}. At the light of corollary II. 6[ this means that a 
combinatorial study of the matrix entries of the operator R + L + S appearing in the 
work of Sukumar & Hodges (2007) is easiest when (a — l)/2 or (a + l)/2 belongs to 
{0, 1}. As a consequence of (11. 6p we obtain 

(eo, {L + S + R)^eo) = ( for a G (-1, 1]. (3.5) 
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The case a = —1 may be handled using remark [TTTl Similarly, for ei we have 

(ei,(L + ^ + i?)"ei) = (^^^^ for all a e [-1,1]. (3.6) 

Another combinatorially interesting case is c = where proposition l3.2| gives p[k, 0]„(x) = 
Mn{x; 0, 2k). As a consequence, the operator L + R appearing in the work of Sukumar 
& Hodges (2007) satisfies 

(eo,(i: + i?)2"eo) =;x2„(0,(« + l)/2) for «G (-1,1]. (3.7) 

The case a = —1 may again be handled using remark [TTTl Similarly, for ei we have 

(ei, (L + i?)2"ei) = /i2n(0, (a + 3)/2) for all a G [-1, 1]. (3.8) 

Here fi2n{^,v) stands for the (2?T,)th moment associated to the Meixner polynomials of 
the second kind with parameters (0,//). For 1] = 1, fll.lOp takes a very simple form, as 
this is noted in Example 21 of ch. II by Viennot (1983). Thus, by applying corollary 13. II 
to c = and = 1/2, we obtain a combinatorial explanation for the first half of (11.221) . 

Remark 3.4. Viennot's (1983) work offers a nice model for the moments /i2n(0, 1) by 
further restricting the appropriate "Laguerre stories". This corresponds to enumer- 
ating certain zig-zag permutations. For other values of rj we could not avoid avoid 
considering "weighted stories" in Viennot's (1983) model. The case a = in the work 
Sukumar & Hodges (2007) does not correspond to r] = 1, yet Sukumar & Hodges (2007) 
find a simple, non-weighted model by considering the enumeration of transposition zig- 
zag classes. It is an interesting question for future research whether this enumeration 
question may be related to Viennot's (1983) theory by constructing a highly nontriv- 
ial bijection, similarly to Viennot's (1983) bijection between all permutations and his 
"Laguerre stories" . 



For other values of rj we may use lemma 11.31 which gives 

j.n 

(eo,vr+(5-C)'"eo)=i5;S^ where sec^^t) = J] i^f ) - (3.9) 

n>0 

In particular, substituting k = 1 into (13.91) yields the second half of (11.221) . Indeed, as 

/r)\ 

noted in eq. (3.2) of Sukumar & Hodges (2007), we have E^n ' = Ta^+i. Using for 
an arbitrary parameter k allows us to rewrite (13. 7p and (13. 8p as 

(eo, {L + Rf'eo) = E^^^^^'^'^^ for a G (-1, 1] and (3.10) 

(ei, (L + i?)2"ei) = E'i^^^^'^^ for all a G [-1, 1]. (3.11) 

Finally, in the case |c| > 1, leading to Meixner polynomials of the first kind, we 
may use ([H]) to find (eo, (vr^(-C) + c • T^t{H) + 7r+(5))'"eo). We omit the details for 
this case and mention only one, combinatorially interesting, subcase. For k = 1/2 and 
c = 3/-\/8 the parameter (3 of the associated Meixner polynomial becomes 1, whereas 
the parameter "c" becomes 1/2. As it is explained in ch. II, §8 of Viennot's work 
(1983), the unitary Meixner polynomials of the first kind {m„(a;; 1, l/2)}„>o are the 
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Kreweras polynomials whose moments /i„ are the numbers of bicolored permutations of 
{1,2, ... ,n} where the two colors are used to color the descents of each permutation. 

What we have seen thus far is that the application of corollary 12.61 to positive dis- 
crete series representations of 1) is linked to the use of Laguerre polynomials and 
of Meixner polynomials of both kinds. The same orthogonal polynomial sequences need 
to be used to apply theorem 11.21 to the same representations. The only difference is 
that, in each case, instead of the moments we would need to use the inverse polynomi- 
als of the same OPS. We may find these inverse polynomials using fll.lSp and (11.141) . 
For brevity's sake, in the rest of the section we only outline how to find the inverse 
polynomials of Laguerre and Meixner polynomials, and leave the adaptation of these 
data to theorem 11.21 to the reader. 

The Laguerre polynomials {L^\x)}n>o satisfy (29) in ch. II of Viennot's (1983) work: 

fn / f \ 

J2^n\^)- = (1 + t)-"-^exp , (3.12) 

n>0 ^' V + / 

thus we must set f(t) = (1 + t)^""^^ and g(t) = t/(l + t) in (11.131) . As indicated in 
table 4 of ch. Ill by Viennot (1983), this implies g^^^^{t) = t/{l — t). Thus the inverse 
polynomials of {L^n\x)}n>o are given by 

r / 1 ft 



Hence we have 



5^Q„(a:)-=U-^ exp x^). (3.13) 

n>0 ■ \ / \ 



mw-||:(-i)'(-7^)(„_l,>--"-- 



{-ly-'^nl f-d-a-1 



d\ V n — d 



implying 



/^n,d = ( ^ ) (a + 1 + C?)n-d- (3.14) 



The generating function for the unitary Meixner polynomials {m„(x; /?, c)}„>o may 
be found in table 4 of ch. Ill by Viennot (1983): 

n.>0 ^ ^ ^ ^ 

thus we must set f{t) = (1 + tc/(l - c))~'^ and g{t) = ln(l - t + c) - ln(l - c + ct) 
in (11.131) . As indicated in table 4 of ch. Ill by Viennot (1983), this implies g^~^^{t) = 
(1 — c)(exp(t) — 1)/(1 — cexp(t)), and the inverse polynomials of {m„(x; /5, c)}„>o are 
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given by 



n>0 



= (-^^) " exp f ^(l-c)(expW -i: I , (3,15) 

^ ^n! Vl-cexp(t)y l-cexp(t) ' ^ ^ 



The Meixner polynomials of the second kind {M„(x; ,5, ?7)}„>o satisfy (69) in ch. II 
of Viennot's (1983) work: 

^M„(x;5,77)^ = ((1 + (5t)2 + t2)-'7/2gxp fxarctan (^^) ) , 

n>0 ^' ^ V + / / 

thus we must set f{t) = ((1 + 5^)^ + t2)-'?/2 ^nd ^(t) = arctan(t/(l + 6t)) in ffTTBD . 
As indicated in table 4 of ch. Ill by Viennot (1983), this implies g^~^\t) = tan(t)/(l — 
5tant), and the inverse polynomials of {Mn{x] , S,r])}n>o are given by 

v:^^ , ,r / sec(t)2 V^' f tan(t) \ , , 

In particular, we have 

;x„_^ = _[r]sec(t)nan(t)'^ for 5 = 0. (3.17) 

(Jj. 

4. Eigenvectors of closable tridiagonal operators 

In this section we return to general tridiagonal operators on the complex Hilbert space 
£^(Z>o), but it will be convenient to assume that they are closable. This assumption 
was not necessary in section O where we were only interested in calculating the matrix 
entries in the first column for the powers of tridiagonal operators. Formally, we are still 
able to state theorem 14.11 below without assuming closability. However, the theorem 
implies that the eigenvectors of tridiagonal operators of the stated form must arise as 
infinite linear combinations of the basis vectors. In all applications, it seems easier 
to show that such vectors belong to the domain of the operator when the operator is 
closable. 

As noted at the end of section 2.1 by Groenevelt & Koelink (2002), for any (positive 
or negative) discrete series representation tt of su{l, 1), the operators vr(i?), vr(C) and 
7t{H) are closable. As a consequence, the operators operators Tc defined in (13.11) are 
closable. They are even self-adjoint, i.e., Jacobi operators. We refer the reader to 
Koelink (2004) for the spectral theory of Jacobi operators. This theory was used by 
Van Der Jeugt (1997), Koelink & Van Der Jeugt (1998) and Groenevelt & Koelink 
(2002) to obtain a complete description of the spectra of the operators Tc. The main 
theorem of this section generalizes only the necessity part of their results. Even this 
generalization suffices to highlight the fact: the (potential) eigenvectors of any (closable) 
tridiagonal operator T may be expressed using the same OPS, whose moments and and 
inverse polynomial coefficients may be used to express the first columns in the matrices 
of powers of T. 
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Theorem 4.1. Let T be a tridiagonal operator with LDU -decomposition T = L + D + U 
such that InUn 7^ holds for all n > and let z ^ C any complex number. If z is an 
eigenvalue ofT then the associated eigenspace is one dimensional, generated by 

_ \^ Pn{z) 

„>0 lli=0 

Here {pn(a^)}n>o is the sequence of monic polynomials defined by ( li. il) and ( (i.^) . sat- 
isfying bn = dn and A„ = Un-lln-l- 



Proof. A vector of the form X]n>o ^" ' eigenvector of T associated to the eigenvalue 
z only if 

d^hQ + Uohi = zHq and (4.1) 

^n-l^n-l + dnhn + Unhn+l = zkn holds for U > 1. (4.2) 

Introducing po := Hq and 

Pn := hn ■ UqUi ■ ■ ■ for n > 1 

we may rewrite (14.10 as 

Pi = {z- do)po, (4.3) 
while multiplying both sides of (14. 2 p by UqUi ■ ■ ■ Un~i yields the equivalent equation 

Pn+l = {Z- dn)Pn " for U > I. (4.4) 

If Po = then (14.31) and (14.41) implies Pn = for all n. We obtain a nonzero eigenvector 
only if po 7^ and then, without loss of generality we may assume po = 1. Furthermore, 
(14.41) implies that p„ must be obtained by substituting z into a monic polynomial 
sequence {pn{.x)}n>o satisfying Favard's recursion formula with 6„ = d^ and A„ = 
Un-iln~i for n > 1. Therefore any eigenvector associated to z is of the form stated in 
the theorem. □ 

Remark 4.2. In the special case when the tridiagonal operator T is a Jacobi operator, 
T is self-adjoint and thus has only real eigenvalues. 



The connecting coefficients {A„}„>i resp. {bn}n>o for the polynomials {p„(a;)}„>o used 
in theorem 14.11 above are the same as the ones used in corollary 12.61 and theorem 12. 8[ 
Thus, the same OPS may be used to calculate the first column of matrix entries of the 
powers of T as in the calculation of its eigenvectors. In particular, for any operator Tc 
defined in (13.11) we must use the appropriate OPS {p[k, c]n(a;)}„>o from proposition 13. 2[ 
This leads to the use of Laguerre polynomials, and Meixner polynomials of both kinds, 
as in the above cited results of Van Der Jeugt (1997), Koelink & Van Der Jeugt (1998) 
and Groenevelt & Koelink (2002). 

We conclude this section with an observation regarding the case c = and k = 1/2, 
leading to p[l/2,0]{x) = M„(x;0,l), by proposition 13.21 A combinatorial model for 
the coefficients of the powers of x in M„(x;0, 1) is mentioned in Example 21 of ch. 
II by Viennot (1983). Up to sign, the coefficient of in M„(x;0, 1) is the number 
of permutations of {l,2,...,n} having j odd cycles. Thus, for this specific choice 
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of parameters, not only the selected matrix entries of the powers of Tq, but also the 
description of the eigenvectors has a combinatorial interpretation. Recall that the same 
selection of parameters leads to the first part of (11.221) . which is equivalent to the first 
part of (11.211) that also caught the attention of Hodges & Sukumar (2007). Therefore, 
if someone wanted to build a "purely combinatorial model" of a quantum oscillator in 
the future, starting with a = 1 in the model proposed Hodges & Sukumar (2007) or 
with /c = 1/2 in the model proposed by Klimyk (2006), and focusing on understanding 
the action of the moment operator seems a very reasonable first step. 



Up to this point, our apphcations of theorems 12. 5[ 12.8 1 and 14.11 involved self-adjoint 
operators from discrete series representations of su{l, 1). In this final section we outline 
two examples indicating that that these theorems may also be useful in other applica- 
tions where the operator is not self-adjoint or where it belongs to the representation of 
an algebra that is is not sm(1, 1). 

Our first example is tt'^{B + C) where tt^ is a positive discrete series representation of 
su{l, 1) for some k > 0. In analogy to corollary 13. 11 the matrix entry (cq, 'n'^{B + C)"^eo) 
is the mth moment of the functional associated to the OPS {fn''\x)}n>o given by (II. ip 
and (11.21) where 

bn = and \n = —n{n -|- 2A; — 1). 
It is easy to verify that we have 



where the polynomials M„(x; 0, 2k) are Meixner polynomials of the second kind. Thus, 
for k G |Z>o, the polynomials fn^\x) are identical to the polynomials studied by 
Carlitz (1959). The applications of theorems 12.51 and 12.81 are no different from the self- 
adjoint setting. By theorem 14. H a complex number z G C is an eigenvalue 7r^(-B + C) 
only if the expression 



represents a vector belonging to £^(Z>o). The the discussion of the convergence issues 
is analogous to the case of vr^(i? — C) that was worked out by Koelink & Van Der Jeugt 
(1998): here we obtain that all eigenvalues belong to iM. This is not surprising, since 
m^i^B + C) is self-adjoint. However, the way were able to handle the example indicates 
that the validity of our main results is not restricted to self-adjoint operators. 



Our second example is the sum A -|- A''' where A and are the operators defined in 
eq. (2.1) of Sukumar & Hodges (2007): 



5. Beyond self-adjoint operators and beyond sm(1, 1) 



f^^\x) = {-irMr,{ix-{),2k) 



(5.1) 




Ae2n 



!n 





Ae2n+1 



y/2n + 1 + a ■ 62; 
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The operator A + is self-adjoint, but it belongs to a representation of a Lie superal- 
gebra properly containing su{l, 1). Theorems 12.51 [2781 and I4.ll are obviously applicable, 
and the OPS {pn{x)}n>o appearing in all of them is given by (11.11) and (11. 2p where, 

bn = 0, A2ri,+i = 2n + 1 + a, and = 2?2. 

We obtain the generalized Hermite polynomials, see eq. (2.46) of ch. V in Chihara (1978). 
Choosing a = yields the usual Hermite polynomials. Sukumar & Hodges (2007) note 
that this choice corresponds to the standard harmonic oscillator. For nonzero a, the 
parity-dependent nature of the quantum algebra defined by Sukumar & Hodges (2007) 
is reflected in the parity-dependent nature of the recursion defining the generalized 
Hermite polynomials. 

Acknowledgements 

I wish to thank Mireille Bousquet-Melou for acquainting me with Viennot's (1983) 
work and two anonymous referees whose suggestions greatly improved the content and 
presentation of this paper. Many thanks to Erik Koelink and Alan Lambert who 
promptly answered many of my questions. This work was supported by the NSA grant 
^ H98230-07-1-0073 and essentially completed while the author was on reassignment 
of duties sponsored by the University of North Carolina at Charlotte. 

References 

[I] Carlitz, L. 1959 Bernoulli and Euler numbers and orthogonal polynomials. Duke Math. J. 26, 
1-16. 

[2] Chihara, T. S. 1978 An introduction to orthogonal polynomials. New York-London-Paris: Gordon 
and Breach Science Publishers. 

[3] Groenevelt, W. 2004 Tensor product representations and special functions, doctoral thesis, Delft 
University Press. 

[4] Groenevelt, W. & Koelink, E. 2002 Meixner functions and polynomials related to Lie algebra 
representations, J. Phys. A 35, 65-85. 

[5] Hodges, A. & Sukumar, C. V. 2007 Bernoulli, Euler, permutations and quantum algebras. Proc. 

R. Soc. A 463, 2401-2414. 
[6] Klimyk, A. U. 2006 The su(l,l)-models of quantum oscillator. Ukr. J. Phys. 51, 1019-1027. 
[7] Koekoek, R. & Swarttouw, R. F. 1998 The Askey-scheme of hypergeometric orthogonal polynomials 

and its (/-analogue, http://fa.its.tudelft.nl/~koekoek/askey/ 

[8] Koelink, E. 2004 Spectral theory and special functions, Laredo lectures on orthogonal polynomials 
and special functions, 45-84, Adv. Theory Spec. Funct. Orthogonal Polynomials, Hauppage, NY: Nova 
Sci. Publ. 

[9] Koelink, E. & Van Der Jeugt, J. 1998, Convolutions for orthogonal polynomials from Lie and 
quantum algebra representations, SIAM J. Math. Anal. 29, 794-822. 

[10] Stanton, D. 2003, Enumeration and special functions, in "Orthogonal polynomials and special 
functions (Leuven 2002)", 137-166, Lecture Notes in Math. 1817, Berlin: Springer. 

[II] Sukumar, C. V. & Hodges, A. 2007 Quantum algebras and parity-dependent spectra, Proc. R. 
Soc. A 463, 2415-2427. 

[12] Van Der Jeugt, J. 1997 Coupling coefhcients for Lie algebra representations and addition formulas 
for special functions, J. Math. Phys. 38, 2728-2740. 



20 



G. HETYEI 



[13] Viennot, X. 1983 Une theorie combinatoire des polynomes orthogonaux generaux. Montreal: Lec- 
ture Notes LACIM, UQAM, 219p. Publication du LACIM, Universite du Quebec a Montreal, 1984, 
reed. 1991. Available online at the author's personal website |http : //web . mac . com/xgviennot| 

[14] Vilenkin, N. Ja. & Klimyk, A. U. 1991 Representation of Lie groups and special functions, vol. 1. 
Dordrecht: Kluwer. 

E-mail address: ghetyei@uncc.edu 

Department of Mathematics and Statistics, University of North Carolina at Char- 
lotte, Charlotte, NC 28223, USA 



